Search for: All records

This study introduces novel constraints on the free streaming of thermal relic warm dark matter (WDM) from Lyman- $\alpha$forest flux power spectra. Our analysis utilizes a high resolution, high redshift sample of quasar spectra observed using the HIRES and UVES spectrographs ( $z=4.2–5.0$). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters including the free streaming of dark matter, cosmological parameters, the thermal history of the intergalactic medium, and inhomogeneous reionization to establish lower limits on the mass of a thermal relic WDM particle of 5.7 keV (at 95% CL). This result surpasses previous limits from the Lyman- $\alpha$forest through reduction of the measured uncertainties due to a larger statistical sample and by measuring clustering to smaller scales ( $k_{\max }=0.2{\mathrm{km}}^{-1}\mathrm{s}$). The approximately two-fold improvement due to the expanded statistical sample suggests that the effectiveness of Lyman- $\alpha$forest constraints on WDM models at high redshifts are limited by the availability of high quality quasar spectra. Restricting the analysis to comparable scales and thermal history priors as in prior studies ( $k_{\max }<0.1{\mathrm{km}}^{-1}\mathrm{s}$) lowers the bound on the WDM mass to 4.1 keV. As the precision of the measurements increases, it becomes crucial to examine the instrumental and modeling systematics. On the modeling front, we argue that the impact of the thermal history uncertainty on the WDM particle mass constraint has diminished due to improved independent observations. At the smallest scales, the primary source of modeling systematic arises from the structure in the peculiar velocity of the intergalactic medium and inhomogeneous reionization. Published by the American Physical Society2024 

ABSTRACT We present the Sherwood–Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations, and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-α forest. Our main findings are: (i) consistent with previous studies patchy reionization causes large-scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-α forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-α forest, such as flat-bottom or double-dip absorption profiles. 

ABSTRACT We measure the mean free path ($$\lambda _{\rm mfp,H\, \small {I}}$$), photoionization rate ($$\langle \Gamma _{\rm H\, \small {I}} \rangle$$), and neutral fraction ($$\langle f_{\rm H\, \small {I}} \rangle$$) of hydrogen in 12 redshift bins at 4.85 < z < 6.05 from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modelled by post-processing simulations from the Sherwood suite using our new code ‘EXtended reionization based on the Code for Ionization and Temperature Evolution’ (ex-cite). ex-cite uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $$\Gamma _{\rm H\, \small {I}}$$ fluctuation models. Our simulation with ex-cite shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within 1σ uncertainty. We measure the three parameters by forward-modelling the  Lyα forest and comparing the effective optical depth ($$\tau _{\rm eff, H\, \small {I}}$$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modelling parameters, observational systematics, and cosmic variance. Our best-fitting parameters show significant evolution with redshift such that $$\lambda _{\rm mfp,H\, \small {I}}$$ and $$\langle f_{\rm H\, \small {I}} \rangle$$ decreases and increases by a factor ∼6 and ∼104, respectively from z ∼ 5 to z ∼ 6. By comparing our $$\lambda _{\rm mfp,H\, \small {I}}$$, $$\langle \Gamma _{\rm H\, \small {I}} \rangle$$ and $$\langle f_{\rm H\, \small {I}} \rangle$$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best-fitting parameter evolution is consistent with a model in which reionization completes by z ∼ 5.2. Our best-fitting model that matches the $$\tau _{\rm eff, H\, \small {I}}$$ distribution also reproduces the dark gap length distribution and transmission spike height distribution suggesting robustness and accuracy of our measured parameters.